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Abstract 


It is shown that , by means of canonical operator approach, the Ward-Takahashi identity (WTI) 
at finite temperature T and finite chemical potential jj, for complete vectorial vertex and complete 
fermion propagator can be simply proven, rigorously for Quantum Electrodynamics (QED) and ap¬ 
proximately for Quantum Chromodynamics (QCD) where the ghost effect in the fermion sector is 
neglected. The WTI shown in the real-time thermal matrix form will give definite thermal constraints 
on the imaginary part of inverse complete Eeynman propagator including self-energy for fermion and 
will play important role in relevant physical processes. When the above inverse propagator is assumed 
to be real, the thermal WTI will essentially be reduced to its form at T — fj, = 0 thus one can use it 
in the latter’s form. At this point, a practical example is indicated. 
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Ward-Takahashi identity (WTI) is a relation of fundamental importance in field theory, including 

in a thermal gauge theory. For instance, in thermal QCD, when discussing the Schwinger-Dyson equation 
for guark self-energy and the effective potential for quark propagator, one needs a non-perturbative WTI 
for complete vectorial vertex and complete fermion propagator, since in these cases the used fermion 
propagators are complete ones without perturbative corrections. If we use the real-time formalism of 
thermal field theory, then the WTI should be expressed in a thermal matrix form. On the WTI in the 
real-time thermal field theory, many researches have been made, e.g. a general demonstration of WTI 
from canonical operator approach of fields |S], a formal discussion of WTI from the path integration 
approach ^ and an explicit check of WTI from a perturbation theory etc.. However, these researches 
did not involve the explicit matrix expression of the above non-perturbative WTI. Hence, for practical 
purpose, it is interesting to make a deep-going investigation of the non-perturbative WTI. The rigorous 
thermal WTI in a gauge theory, in principle, may be derived from Becchi-Rouet-Stora (BRS) invariance 
1^. However, when the problem is confined to the WTI for vectorial vertex in a vectorial gauge theory, 
we can alternatively show a more simple approach to derive it. In this paper, we will indicate that, 
following the same approach to obtain the WTI in zero temperature QED, the non-perturbative WTI at 
a finite temperature T and a finite chemical potential fi can be proven simply by means of the operator 
formalism of the fields, rigorously for QED and approximately for QCD. Then the feature of the resulting 
WTI and its relation to practical application will be analyzed further. 

First let us consider the case of temperature QED. In QED, gauge invariance amounts to electroma- 
netic current conservation, hence the WTI at finite T and ft can also be derived from the latter. In the 
real-time formalism of thermal field theory [J, the conserved Noether vectorial currents j^\x) may be 


expressed by 


j\^\x) = 'ip^'^\x)-f^'ip^’'\x), r = 1,2 


( 1 ) 


which obeys 


9"jT(^) = 0, r=I,2, 


( 2 ) 
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where are the fermion field operators, r = 1 and 2 respectively denote physical and thermal ghost 

fields (here the attribute ’’thermal” has been added so as to distinguish the ghost in thermal field theory 
from the ghost in a gauge theory). Now define 


jv{x) = X! 

r^l,2 


(3) 


then it is also conserved, i.e. 

d'' ju{x) = 0. 

We first prove that conservation of the current (x) will lead to the equality 

d^Tj^{x)'tlj^^\xi)ip^‘'\x2) = [-S'^{x - Xi) + 6 '^{x - X2)]T'ip^''\xi)jp^‘'\x2),r,s = 1,2, 


(4) 

(5) 


where represents only the partial derivative for x and T is the time-order operation. In fact, by means 
of the definition of the T product and Eq.(4), we can obtain 

d:ij„{x)^p^'^\xi)^^^\x2) = 5{x°-x\)l[j^{x),^lj^'-\xi)\^^^\x2) 

-b(5(a;° - xl)lilj^'^\xi)[j^{x),'il)^‘'\x 2 )\. 

Considering the fact that the physical fields and the thermal ghost fields commute (or anticommute) each 
other [7], we will always have 


[jo(x),V’(’')(xi)] 

[jo{x),^^^\x2)] 


-S^{x - Xi)ip^'''>(xi), 
S^{x - X2)'tp^^\x2), 


and then are led to Eq. (5) being proven. 

Now taking thermal expectation values of the both sides of Eq. (5) and defining the thermal three-point 
function 

G"(x;xi,a;2) = {0\Tj^{x)'ip^^\xi)ip^''\x2)\0)T (6) 

and the thermal propagator 

G^%XI,X2) = (0|TV^W(xi)^(^)(cr2)|0}T (7) 

where (0| ■ • ■ |0 )t represents thermal expectation value, then we will get 

d^Gl''{x]Xi,X 2 ) = [-S‘^{x - xi) + 6‘^{x - X 2 )]G'’''{xi - X 2 ), r,s = 1,2. (8) 

In Eq. (8) the translation invariance of G'^‘‘(xi,X 2 ) has been considered. The Eourier transformation of 
G0®(a;; Xi,X 2 ) and G’'®(xi — X 2 ) are denoted respectively by Glf{q;pi,p 2 ) and iS'^‘^{pi) or iS''~‘^{p 2 ), where 
Glf‘{q;pi,p 2 ) is the complete vectorial vertex and iS''"’^(p) is the complete fermion propagator. Then the 
Eourier transformation of Eq.(8) will lead to that 

iq''Gl^q-p,,P 2 )=lS'''^{p 2 )-lS'^^{pi), q = P 2 -Pl. (9) 


On the other hand, G”'’(x;xi,X 2 ) is the three-point function with external propagators (legs), so it can 
be written by 

G(;^(x;xi,X 2) = ^5'”‘(xl - x)Tl*\x)iS'*'^ix - X 2 ). (10) 

The Fourier transformation of Eq.(lO) leads to that 

Glfiq]Pi,P 2 ) = iS'^\pi)Tl^{p 2 ,Pi)iS'*'‘'{p 2 ), q = P 2 - Pi- (H) 

Substituting Eq.(ll) into Eq.(9) we obtain 

tS'^\pi)tq''Tl*'iP2,pi)tS'^'^{p2) = tS'^^p2)-tS'^^{pi). ( 12 ) 

In the form of thermal matrix, Eq. (12) can be expressed by 

iS'{pi)iq''t'^{p 2 ,pi)iS'{P 2 ) = iS'{p 2 ) - iS'ijpi). (13) 


2 




The thermal matrix propagator iS'{p) is defined by |7| 

where 

£r _ f cos 9p sm9p\ 

p sin9p cos dp J 

is the thermal transformation matrix with fermionic chemical potential /i, /3 = 1/T , 

P e0(p°-p-) + 1 e^^~P°+P'^ + 1 

and iSp^p) and —iSp{p) are respectively the complete Feynman propagator for fermion and its complex 
conjugate, but the complex conjugate does not include the 7 matrix, in the light of the convention 
in the real-time thermal field theory |^. The thermal vertex matrix T^{p 2 ,pi) is defined by 


i'i^{P2,Pi) 


{^l^{P2,Pl) ^i^{P2,Pi)\ 

\^f{P2,Pl) ^f{P2,Pl))' 


The matrix equation (13) can be transformed to 


= [5'(p2)]-i - [ S ' ip ,)]-\ 


(14) 


Eq.(13) or Eq.(14) is the WTI in thermal matrix form at finite T and /i in QED. They are the rigorous 
expressions of the thermal WTI for complete vectorial vertex and complete fermion propagator and 
obviously, should also be valid to each order of the perturbation theory. 

The WTI at finite T and p expressed by Eq.(13) or Eq.(14) which was derived from the Abelian 1/(1) 
current conservation in QED can not be simply generalized to the non-Abelian SUc{i) case in QCD. 
In QCD the corresponding identity should be the Slavnov-Taylor identity (STI) coming from invariance 
under BRS transformation. However, as seen from the case of STI at T = /r = 0, if one neglects the effect 
of ghosts in the fermion sector inmni, and this amounts to replacing the composite fermion-ghost vertex 
and the ghost propagator in the STI by their bare or the lowest order values, the STI for fermion-fermion- 
gluon vertex will be reduced to the form of the WTI in QED, only a SUc{3) generator matrix must be 
multiplied in the both sided of the identity m The main purpose of this approximation is for simplifying 
practical calculations e.g. making the derived Schwinger-Dyson equation of fermion self-energy become 
tractable nmnnn. We will prove that the approximative form of the STI can be obtained simply by the 
conservation of the SUc{S) vector current, hence neglecting the ghost effect in fermion sector only implies 
that the local QCD gauge invariance now becomes the requirement that the color current is conserved 
at fermion-fermion-gluon vertex. In fact, when the ghost fields are neglected, the Lagrangian of QCD 
including the gauge fixing term is invariant under a global SUc{S) transformation and this will lead to 
conservation of the color current 


a, 6 ,c = I, • • •, 8 . (15) 

where A“(a = 1, • • ■ , 8 ) are the SUc{3) generators, are the structure constants of SUc{3) group, and 
and are the gauge fields and their field strengths. Eq.(15) was derived in the approximation that 
the ghost fields are completely neglected. However, we will use it only in the case of fermion-fermion- 
vertex, i.e. the ghost effect is neglected only in fermion sector. Thus by means of Eq.(15) and the same 
demonstrations as that made in QED, we may obtain the WTI in QCD with the similar form to the one 
in QED. 

Assuming the above approximation is kept at finite T and p, then one will be able to derive the thermal 
WTI for fermion-fermion-gluon vertex in the same way as the above in QED. As a non-Abelian gauge 
theory, QCD is normally described by path integral approach; however, when considered as a field system 
with constraints, it may also be discussed in principle by canonical operator approach. In particular, for 
derivation of the WTI for fermion-fermion-gluon vertex, one mainly uses the communicators between the 
fermion fields and does not involve any communicators between the gauge fields. Hence the operator 
approach is completely feasible in this problem. 
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To derive the thermal WTI for fermion-fermion-gluon vertex in the approximation of neglecting the 
ghost effect, we may replace the Abelian vector current (3) by the non-Abelian SUc{3) vectorial current 


r^l,2 





(16) 


where for simplicity, we will consider only a single flavor of fermions. By means of the current conservation 
d''j^{x) = 0 and the similar derivation leading to Eq.(5) we can prove that 


+S\x-X2)T4"\^i)^^'^''i^2){X‘^)k^, r,s=l,2, (17) 

where i and j denote the color components of the spinors and and we have considered the fact 
that 


Noting that the thermal propagators 


(0|T^;’’^(xi)'i/'^®^'^(a;2)|0)7’ = 0 , when I ^ j, 

(0|Ti/)|’’^(xi)V>^®^^(a;2)|0)r = 0 , when i^k, 

though they may be non-zeros when r ^ s and I = j or i = k. Hence the thermal expectation value of 
Eq.(17) will become 


5 x(0|Tj“(x)'i/'|’'^(xi)'i/'^®^'’(a;2)|0)T = (A“)/[- 5 '‘(x - xi) + 6‘^{x - X2)](0|T'i/'^’'na;i)V’^'*Ha;2)|0)T, r,s= 1 , 2 . 

(18) 

Since the SUc{3) generator matrix A“ has been factorized out in the right-handed side of Eq. (18), we 
can define 

(0|Tj“(a;)f/'-’'^(a;i)'!/;(®)'^(x2)|0)T = (A“)/G;;®(x; Xi, a;2). (19) 

We indicate that since the gauge field sector of j“(x) has no contribution to the discussed three-point 
Green function, the G((®(x;xi,X2) defined by Eq.(19) is actually the fermion-fermion-gluon vertex. As a 
result of Eqs.(18)-(19) and Eq.(7), we will have 

d^{X°‘)/Gl%x-,xi,X 2 ) = {X°-)X[-6'^{x - xi) + S*{x - X2)]G’'®(xi - X2), r, s = 1,2. (20) 


Eliminating A“ from the both sides of Eq.(20) we will come to Eq.(8) once again. This implies that, for 
the fermion of single flavor and single color, when the ghost effect is neglected, the WTI for fermion- 
fermion-gluon vertex in QCD is identical to that in QED, i.e. the thermal WTI (13) or (14) are also 
valid in QCD in the approximation that the ghost effect in fermion sector is neglected. It is emphasized 
that the thermal WTI (13) or (14) rigorously for QED and approximately for QCD are proven by the 
operator approach here and it seems that the similar proof did not appear in the literature before. 

By means of the explicit matrix expression of Eq.(14) 


9'T“(p2,pi) g''ri2(p2,pi)\ 



q''Tf{p 2 ,Pi) q''Tf(j) 2 ,pi) j 

Sp^{p) - sin'^ 0p[Sp^{p) - S*p, \p)], 

cos0pSm0p[Sp^{p) - S*p ^{p)] \ 


_g-/3,i/2 cos6>pSin6»p[5'^^(p) - S'^“\p)] 

1 9 1 , - (P2 - 

-^Pi), 

-S*p-\p) - sin^ 0p[Sp\p) - ^r'b)] /p=p. 

(21) 


we may obtain respectively from its (11), (22) and (12) or (21) components the equations 
q''rl\p2,pi) = sp\p2) - sp\pi) - sin^ 0pASp\p2) - S*p-\P2)] + Sin2 0^^ [Sp\pi) - Sp-\pi)], (22) 

r^^(P2,Pi) = -[r“(P2,Pi)]* (23) 
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and 


q''K^iP2,Pi) = -el^^q''Tl^{p2,pi) = cos9p^siii9p^[Sp^(j)2) - S*p ^(^ 2 )] - (P2 ^ Pi)- (24) 
To tree-diagram order, we will have 


r“(P2,Pi) = -TI^{p2,Pi) = lu, 
Sf^{p) = S*p~'^{p) - m, 

and the thermal WTI (21) becomes 

/ q^lu 0 \ ^ ( ^ 2 - A ® ^ 

0 -q''-l^ / \ 0 -(^ 2 - A)) 


(25) 


which is obviously valid. In the limit of T = p = 0, we only need to consider the (11) component of 
Eq.(21). Since 

sin0p = 0, when T = /i = 0, 


we obtain that 

g‘'r“(p2,Pi) = Sp^(j)2) - Sp^ipi) 


(26) 


which, with Tl}(p 2 ,pi) being identified with the physical vertex, exactly reproduces the standard form 
of the WTI at T = /r = 0. 

In the general case of finite T and /i, we see that Eq.(22) contains the thermal correction terms, 
i.e. the last two terms in its right-handed side, and Eq.(24) represents definite thermal constraints on 
T]^{P 2 ,Pi) and r^^(p 2 ,Pi)- These corrections and constraints exist only at finite T and p and always 
related to a factor [5')J^^(p) — S'^~^(p)], i.e. the imaginary part of Sp^{p). Since iSpip) is the complete 
Feynman propagator for fermion, so Sp^{p) will contain complete fermion self-energy and could have an 
imaginary part. If we write 


iSpip) A(p^) ^ — B{p'^) -I- ie 


(27) 


then the factor 

Sp\p) - S*p-\p) = [A{p^) - A*{p^)] ji- [B{p^) - B*{p^)] (28) 

will be related to only the imaginary parts of the functions A{p'^) and B{p'^). Therefore, Eqs.(22)-(24) will 
be thermal constraints on the imaginary part of Sp^{p). When one researches into the physical processes 
involving the imaginary part of Sp^{p) or say, fermion self-energy, e.g. discussing the decay rate of 
particle and Landau damping etc. ncni, such constraints would certainly be quite important and must 
be considered carefully, even the calculations are conducted only up to a given order of perturbation 
theory. On the other hand, for some physical processes which do not relate to the imaginary part of 
Sp^{p) e.g. dynamical chiral symmetry breaking and restoration, one only needs to deal with a real mass 
function of fermion which can come from B(p'^)/A{p^). In this case, we can assume that Sp^{p) is a real 
function, i.e. 

Sp^ip) - S*p~^{p) = 0, for any p (29) 


or by Eq.(28), equivalently 


A{p^)=A*{p\B{p^) = B*{p^). 


(30) 


As a result we will have 


r“(P2,Pi) = -[Tf{p2,pi)]* = -Tf{p2,pi) = T^(p2,Pi) 
and the thermal WTI (21) becomes 

q'"T^{p2,pi) 0 \ ^ f Sp\p2) - Sp\pi), 0 

0 -q''T^{p2,pi) J V 0 -[S'^^(p2) - ^^^(Pi)] 

which is essentially identical to the WTI (26) a,t T = p = 0. However, we emphasize that only if Sp^{p) 
is a real function, the WTI at finite T and p is just reduced to its form a,t T = p = 0. 

The above discussions indicate that, when one researches the Schwinger-Dyson equation of fermion 
self-energy and the effective potential for complete fermion propagator at finite T and p in the real-time 
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formalism of thermal field theory, one will be able to use the thermal WTI in its form at T = /i = 0. For 
example, for the 1-type vertex of physical fields, the corresponding thermal WTI can be written from 
Eq.(31) and expressed by 

(P2 -PiY^v{P2.Pi) = A{pI) ^2 - B{pI) - A{p\) -b B{pI). (32) 


To satisfy Eq.(32), a possible r,y(p 2 ,Pi) is 


TYp2,Pi) =ly+ _p^)2 {[^(P2) - 1] ^2 - [A{pI) - 1] ^1 - B{pI) + B{pI)), 


(33) 


where an unspecified transverse pieces have assumedly been neglected. Such selection given by Eq.(33) 
could make practical calculations be simplified greatly if the Landau gauge is taken. In fact, in this 
gauge, the factors in thermal gauge held propagator which are coupled to T^{p 2 ,pi) will have the forms 


D‘'YP2-Pi)r^A‘'^{p2-pi)- 


—l 


A^^(P2 - Pi) = (34) 


(P2-Pi)2 + Ze’ " {p2-pi)^+ie 

and corresponding complex conjugate Zl^^*(p 2 — Pi) and A^^*(p 2 — Pi)- Thus we obtain that 


r,.(p 2 ,pi) 


A"^(P 2 -Pi) 

A"''*(P 2 -Pi) 


A"''(P2-Pi) 
A"^*(P2-Pi) ’ 


(35) 


i.e. effectively only a tree vertex is left. Physically this is a well-known result: in the Landau gauge, 
the gauge held propagator is transverse and the effective coupling will be independent of the longitudinal 
sector of ri,(p 2 ,Pi)- 

In conclusion, it has been shown that, based on vectorial current conservation and canonical operator 
approach of helds, we can give a simple proof of the Ward-Takahashi identity at hnite temperature and 
hnite chemical potential for complete vectorial vertex and complete fermion propagator which is rigorous 
in QED and approximate in QCD where the ghost effect in fermion sector has been neglected. The results 
are expressed in thermal matrix form in the real-time thermal held theory. All the thermal corrections 
to the WTI only relate to the imaginary part of the inverse complete Feynman propagator for fermion 
and will play important role in the physical processes associated with the imaginary part of fermion 
self-energy. When the inverse fermion’s Feynman propagator is assumed to be real, the thermal WTI 
will essentially be reduced to its form at zero temperature and zero chemical potential thus one can use 
it in the latter’s form. At this point, a practical example is also mentioned. 
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